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Abstract 

The physical origin of the nuclear symmetry energy is studied within the relativistic mean field 
(RMF) theory. Based on the nuclear binding energies calculated with and without mean isovector 
potential for several isobaric chains we confirm earlier Skyrme-Hartree-Fock result that the nuclear 
symmetry energy strength depends on the mean level spacing e(A) and an effective mean isovector 
potential strength k(A). A detailed analysis of the isospin dependence of the two components 
contributing to the nuclear symmetry energy reveals a quadratic dependence due to the mean- 
isoscalar potential, ~ e T 2 , and, completely unexpectedly, the presence of a strong linear component 
~kT(T + 1 + e/k) in the isovector potential. The latter generates a nuclear symmetry energy in 
RMF theory that is proportional to E sym ~ T(T+ 1) at variance to the non-relativistic calculation. 
The origin of the linear term in RMF theory needs to be further explored. 
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One of the most important topics in current nuclear physics is to search for the existence 
limit of atomic nuclei, i.e., to determine the nuclear drip line. In this respect, the role 
of the continuum in loosely bound nuclei and, in particular, its impact on the treatment 
of pairing correlations has been discussed to great extent in recent time. However, the 
proper understanding and correct reproduction of the nuclear symmetry energy (NSE) may 
have even greater bearing for masses of loosely bound nuclei and certainly is a key issue in 
the study of exotic nuclei. The very fundamental questions in this respect concern both the 
understanding of the microscopic origin of the NSE strength as well as its isospin dependence. 
The latter issue has attracted recently great attention also in N ~ Z nuclei, see Ref. [1] and 
references therein. 

The NSE is conventionally parametrized as: 

£ sym = a sym {A)T{T + A) (1) 

where T = \T Z \ = \N — Z\/2. The strength of the NSE admits typically volume and surface 
components a sym (A) = a v /A — a s /A 4 ^ 3 and its physical origin is traditionally explained 
in terms of the kinetic energy and mean isovector potential (interaction) contributions i.e. 
a S ym(A) = akin (A) + ai n t (A), respectively [2]. The linear term is found to be strongly model 
dependent and there is a common belief that mean-field models yield essentially only a 
quadratic term A 0. On the other hand, the nuclear shell-model [3-5] or models restoring 
isospin symmetry [6] suggest that A ~ 1. No consensus is reached so far concerning the value 
of A although there is certain preference for A ~ 1. Indeed, experimental masses of nuclei 
with small values of T supports the existence of the linear term [7]. Similar conclusions 
were reached by Janecke et al. [8] based on the analysis of experimental binding energies for 
A < 80 nuclei. One of the most accurate mass formula, the so called FRDM [9] employs a 
value of A 1 but inconsistently admits only a volume-like linear term. Assuming T(T + 1) 
dependence Duflo and Zuker have performed a global fit to nuclear masses obtaining [10] 

. „ 134.4 203.6 p ~\ 7"i 
a sym (A) = — - [MeV]. (2) 

A different view on the origin of the NSE was presented recently by Satula and Wyss. In 
Refs. [11-13] it was demonstrated using the Skyrme-Hartree-Fock (SHF) model that the NSE 
can be directly associated with the mean level spacing e(A) and mean isovector potential, 
E sym = \e(A)T 2 + -k(A)T(T + 1) [11-13]. Surprisingly, the self-consistent calculations 



revealed that the complicated isovector mean potential induced by the Skyrme force is similar 
to that obtained from a simple interaction -k(A)T ■ T, i.e., is very accurately characterized 
by a single strength k(A) [11-13]. This study revealed also that the SHF theory yield in fact 
a (partial) linear term with A ~ k/ (2a sym ) and that this term originates from neutron-proton 
exchange interaction. 

Alongside with the SHF calculation, the relativistic mean field (RMF) theory has been 
used for a large variety of nuclear structure phenomena [14]. Since the RMF theory is based 
on a very different concept from the SHF, it is highly interesting to investigate the structure 
of the NSE in the framework of the RMF theory. 

The details of RMF theory together with its applications can be found in a number of 
review articles, see for example Ref. [15] and references therein, and will not be repeated 
here. The basic ansatz of the RMF theory is a Lagrangian density whereby nucleons are 
described as Dirac particles which interact via the exchange of various mesons [the isoscalar- 
scalar sigma (a), isoscalar- vector omega (u) and isovector- vector rho (p)] and the photon. 
The a and uj mesons provide the attractive and repulsive part of the nucleon-nucleon force, 
respectively. The isospin asymmetry is provided by the isovector p meson. Hence, by 
switching on and off the coupling to the p meson, one can easily separate the role of isoscalar 
and isovector parts of the interaction and study them independently. 

In the nuclei considered here, time reversal symmetry is preserved and the spatial vector 
components of u>, p and A fields vanish. This leaves only the time-like components ou°, p 
and A . Charge conservation guarantees that only the third component of the isovector p 
meson is active. For reason of simplicity, axial symmetry is assumed in the present work. The 
Dirac spinor tpi as well as the meson fields can be expanded in terms of the eigenf unctions 
of a deformed axially symmetric oscillator potential [16] or Woods-Saxon potential [17], and 
the solution of the problem is transformed into a diagonalization of a Hermitian matrix. 

The RMF calculations are performed for the A=40, 48, 56, 88, 100, 120, 140, 160, 164, 
and 180 isobars with the effective Lagrangians NL3 [18], TM1 [19], and PK1 [20]. Our choice 
of the parameterizations is somewhat arbitrary. However, the purpose of this work is not 
to make a detailed comparison to the data but rather to investigate specific features of the 
RMF approach pertaining to the isovector channel. These properties are expected to be 
fairly parameterization independent, in particular that these parameterizations reproduce 
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rather well the equation of state for densities p < 0.2 fm~ 3 [21, 22]. 

The Dirac equations are solved by expansion in the harmonic oscillator basis with 14 
oscillator shells for both the fermion fields and boson fields. The oscillator frequency of the 
harmonic oscillator basis is set to fkuo = AlA^ 1 ^ 3 MeV and the deformation of harmonic 
oscillator basis (3q is reasonably chosen to obtain the lowest energy. Generally speaking, the 
RMF calculation reproduce the experimental binding energy to an accuracy less than 1%. 
For the present study we are mainly interested in the NSE emerging in the RMF theory 
due to the strong (particle-hole) interaction. Hence the Coulomb potentials and the pairing 
correlations will be neglected in the following. The full potential in the Dirac equation is 

V tot = V{v) + /3S(r) = g^°(r) + g p r-p°(r) + (3g a a(r). (3) 

It can easily be separated into isovector and isoscalar components, i.e., V to t = Vi S + Vi V , 
where 

V is (r) = gu u (r)+Pg a a(T), 
V iv (r) = g p f-p°(r). 

The binding energy calculated with the full potential Vtot in Eq. (3) is denoted as Et- 
The energy obtained by switching off the isovector potential, V{ V = 0, i.e. by taking in the 
calculation V to t = V iS) is denoted by E T . In order to single out the impact of isoscalar fields 
on the NSE, we use E T to extract the mean level spacing e(A,T z ) along an isobaric chain 

E T (A,T Z ) - E T=0 (A,T Z = 0) = l -e{A,T z )T 2 . (5) 

The values of e(A,T z ) calculated for the A=48, 88, and 160 isobaric chains of nuclei 
from T z = to the vicinity of the drip line are shown in the left panels of Fig. 1. The 
calculations have been performed for three different parameterizations of the effective La- 
grangian including NL3, TM1, and PK1, respectively, which all yield very similar results. 
For small values of T 2 , strong variations in e{A, T z ) are seen, which are associated with shell 
closures. For larger values of T z , e(A,T z ) become less sensitive to the shell structure and 
its value is stabilized, e(A,T z ) e(A). Note, that the calculated values of e(A) are much 
larger then the empirical estimates for the mean level spacing 53 /A < e emp < 66/ A MeV 
[11, 23-25]. However, after rescaling e*(A) = (m* /m)e(A) by the corresponding effective 
masses m*/m = 0.595, 0.634, and 0.6055 for NL3, TM1, and PK1 respectively, the effective 
mean level spacing e*(A) neatly falls within the empirical bounds (shaded areas), as shown 
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FIG. 1: The mean level spacing e (left) and its counterpart (right) scaled by m*/m for A=48 

(upper), 88 (middle) and 160 (lower) isobaric chains calculated using effective Lagrangians NL3, 

TM1, and PK1 as marked in the figure. The shadowed areas correspond to the empirical mean 
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in the corresponding right panels. Let us also note that with increasing A, i.e., from A = 48 
to A = 160, all curves move toward the upper limit of the empirical data, reflecting the 
decreasing role of surface effects with increasing A, similar to the SHF results [11]. 

When comparing to the results of the SHF calculations in Ref. [11], the following two 
important conclusions can be made: (i) even though the values of s(A) from the RMF 
calculation are much larger than those from the SHF, after the effective mass scaling both 
models generate essentially identical results in agreement to the empirical boundaries for the 
mean level spacing; (ii) the results from the RMF calculations clearly confirm the general 
outcome of Ref. [11], that the isoscalar field generates the NSE of the form eT 2 and that 
its strength is indeed governed by the mean level density rather than the kinetic energy, 
see also [13]. Let us stress that the contribution eT 2 can also be derived analytically using 
the simple iso-cranking model [1, 11, 12]. Note further, that the evidence and conclusions 
gathered from the SHF and the RMF calculations are independent of the iso-cranking model. 

After obtaining the average level density e, we now proceed to calculate the average 
effective strength k of the isovector potential. The effective isovector potential strength 
k is obtained from the binding energy difference between the RMF calculations with and 
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FIG. 2: The average effective strength k of the isovector potential for A=48 (left), 88 (middle) 
and 160 (right) isobaric chains calculated using effective Lagrangians NL3, TM1, and PK1 as 
marked in the figure. Upper, middle and lower panels show the values of k obtained assuming 
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E T — E T — -kT , -kT{T + 1), and —kT(T + 1 + e//c), respectively. See text for further details. 



without the isovector potential for the same nucleus. As explained later, three different 
types of T-dependence of the isovector potential are investigated: 



E t -E t = ^kT 2 and ^kT(T+1) and -kT(T + 1 + e/k). 



(6) 



The resulting effective isovector potential strength k for A= 48, 88, and 160 isobaric 
chains in RMF theory are shown in Fig. 2. Similar to the SHF calculation in Ref. [11], the 
complicated isovector potential along an isobaric chain can be characterized by an effective 
isovector potential strength n(A). However, at variance with the SHF [11], the RMF depen- 
dence on T is best reproduced by a dependence like Et — E~t ~ ~^ K T{T + 1 + e/n), rather 
than T 2 or T(T + 1). Apparently, the linear term in RMF is considerably larger than that 
in the SHF calculation, implying that the total NSE in RMF behaves effectively as: 



£ sym « \eT 2 + ^T(T + 1 + e/k) ~\{e + k)T{T + 1). 



(7) 



The important aspect to note is that in the SHF approximation the linear term originates 
predominantly from the Fock exchange in the isovector channel. In the RMF, which is a 
Hartree approximation, one would therefore expect a kT 2 dependence. In contrast, the 
large slope of the isovector potential strength n(A, T z ), fitted using either kT 2 or kT(T + 1) 



dependence reveals the presence of an effective linear term that is large enough to compensate 
the lack of the linear term in the e proportional term, see Fig. 2. Similar tendencies for the 
mean level spacing e and effective isovector potential strength k are present in all the 10 
isobaric chains we calculated, including A=40, 48, 56, 88, 100, 120, 140, 160, 164, and 180. 

At a first glance, the RMF theory is a Hartree approximation (without exchange term) 
and one does not expect a linear term to be present. On the other hand, the RMF as well 
as the SHF are two particular realizations of the density functional (DF) theory. Moreover, 
due to the large meson masses the relativistic forces should be close to zero range forces and 
both approaches are expected to be rather alike since in this limit the exchange term takes 
the same form as the direct term and should be effectively included within the DF. The 
question that arises is why the RMF is capable to generate the linear term in contrast to the 
SHF approach? In this context it is interesting to observe that within the RMF approach 
the isovector mean-potential is generated by the p-meson field. Hence, its properties are 
defined essentially by a single coupling constant g p , see Eq. (4). In this respect the RMF 
seems to be more flexible than the SHF where isoscalar and isovector parts of the Skyrme 
local energy density functional are strongly dependent upon each other through the auxiliary 
Skyrme force (SF) parameters which are fitted to the data. In the process of detreminig 
the SF, one should therefore balance properly and very carefully the isoscalar and isovector 
data used in the fitting procedure. Indeed, our earlier study [11] shows that the SkO [26] 
parameterization of the SF, which has been fitted to neutron-rich nuclei, has a stronger 
linear term than the so called standard parameterizations but not as strong as that in RMF 
theory. 

We proceed further by investigating the mass dependence of the NSE. Since the total 
NSE in RMF theory behaves effectively like T(T + 1) we extract the NSE strength a syui (A) 
from the difference of the binding energies: 

E syia = E T - E T=0 = a syia (A)T(T + 1). (8) 

To avoid the influence of shell structure we chose two nuclei with large T z for each isobaric 
chain A and calculate a sym (A) simply as an arithmetic mean over such a pair of nuclei. In 
order to compare with the empirical data of Ref. [10], we depict in Fig. 3 the product A*a sym 
as a function of A for all the 10 isobaric chains calculated with the effective Lagrangian 
PK1. The dot-dashed line in Fig. 3 is taken from Ref. [8] revealing the shell structure of the 

7 



> 

CD 



E 
>> 

CO 



CO 

* 

< 



150 



120 



T 



T 



>PK1 

Fitted 

Duflo & Zuker 
Janecke 




200 



FIG. 3: The nuclear symmetry energy coefficient a sym * A (filled circles) extracted from E t 



sym 



o sym T(T + 1) for the A=40, 48, 56, 88, 100, 120, 140, 160, 164, and 180 isobaric chains calculated 
using effective Lag raxigian PK1. The solid line represents fit to the calculations, i.e., A * ttsym — 
133.20 - 220. 27 /A 1 / 3 from Eq. (9). The dashed line is from Eq. (2): A * a sym = 134.4 - 203.6/yl 1 / 3 
[10] and the dot-dashed line is from Ref. [8]. For further details, see text. 

symmetry energy coefficient and the dashed line represents a fit to experimental data given 
by Eq. (2) [10]. The maxima at A=40 and A=56 can easily be seen in our calculation, which 
is in good agreement with Ref. [8]. This behavior is easy to understand since T=0 nuclei in 
the A=40, 56, 100 and 164 isobaric chains are double magic nuclei and hence, more bound, 
resulting in an increase of the symmetry energies for nuclei with T > 0. Still, the average of 
the calculated symmetry energy is quite close to the dashed line, which is fitted to masses 
in Ref. [10]. 

Restricting this analysis to volume and surface terms only, a least-square fit to the cal- 
culated points (filled circles in Fig. 3) leads to a smooth curve: 



a 



(RMF) 
sym 



133.20 220.27 



[MeV], 



A A 4 / 3 

shown as solid line, which is very close to the empirical values (dashed line). 



(9) 



In a similar manner, the volume and surface contributions to the mean level spacing e 
and the average effective strength k are determined from the calculations. In the left panel 
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FIG. 4: Volume and surface contributions to the mean level spacing e (left) and effective strength 

130.42 127.83 



k (right). The fitted line is obtained from the same points as in Fig. 3, i.e., e 
144.77 342.08 

and k - 



A 



A 



of Fig. 4, the mean level spacing e (filled circles) extracted from Ex — Et=o = T^T 2 are 
presented together with the smooth curve obtained from least-square fit: 

130.42 127.83 



A 



[MeV]. 



(10) 



In the right panel of Fig. 4, the average effective strength k (filled circles) extracted from 
E T — E T = -kT(T + 1 + e/k) are presented together with the smooth curve obtained by 



least-square fit: 



144.77 342.08 
k = : 7T7z — [MeV]. 



A 



A 4 / 3 



(11) 



It should be noted that we choose the same nuclei for the fit of e, k, and a sym in the 



least-square fitting. 



It was shown before that for relatively large values of T z both e(A,T z ) w e(^4) and 
k(A, T z ) fa k(A) along each isobaric chain. Let us now fix the value of A and investigate fine 
effects in e(A,T z ), k(A,T z ), and a sym (A,T z ) = -(e(A,T z ) + k(A,T z )) versus T z in order to 
study the response of the isovector potential to changes in shell structure which are naturally 
incorporated in the mean level spacing e(A,T z ). Therefore, we present in Fig. 5 the values 
of e(A,T z ), k(A,T z ), and a sym (A,T z ) calculated using T 2 , T(T + 1 + e/k), and T(T + 1) 
isospin dependencies respectively, for the A = 160 isobaric chain using the parametrization 
PK1 of the effective Lagrangian. To avoid the direct connection between k and e entering 
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FIG. 5: The mean level spacing e, the average effective isovector potential strength k calculated 
assuming either T(T+ 1+e/k) or T(T+ 1 + 1.25) dependence, and the symmetry energy coefficient 
o-sym for the A = 160 isobaric chain in RMF theory with effective Lagrangian PK1. The value of 
1.25 is the average for e/k at A = 160. 

the T(T + 1 + e/k) dependence used to extract k(A,T z ) we also show values of k(A,T z ) 
obtained assuming T(T + 1 + 1.25) dependence, where 1.25 represents the average value of 
e/k. Clearly, there are small variations in e(A, T z ) for T z = 12, 14, 16 and 22 due to the shell 
structure. Whilst the same variations are also obtained for k(A,T z ) at the same value of 
T z , the sum a sym (A,T z ) is very smooth. Apparently, changes in e(A,T z ) and k(A,T z ) cancel 
to large extent. We also see that the variations in k(A,T z ) are not due to the e/k term, 
since similar variations are found for the k(A, T z ) curve calculated using the T(T + 1 + 1.25) 
dependence, see also the curves for T 2 and T(T + 1) as shown in the right panel of Fig. 2. 
Apparently, the variations of k(A, T z ) reflect a direct response to changes in e(A, T z ) showing 
that effectively, the isovector potential and isoscalar potential become closely linked. 

In order to understand the origin of the linear term in RMF theory, self-consistent cal- 
culations for the effective Lagrangian PK1 have been performed including the a and u 
meson-fields only. With the nucleon densities thus obtained, the influence of the p meson 
field has been extracted by switching on the isovector potential in a non-self-consistent way. 
The results for the A=48 isobaric chains are labeled as PK1' in Fig. 6. It is interesting to 
note that just by switching on the isovector potential, the relation T(T + 1 + e/ k) is followed 
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FIG. 6: The comparison of the average effective strength k of the isovector potential for the 
A=48 isobaric chain between self-consistent RMF calculation (PK1) and the corresponding non 
self-consistent (PK1'). See the text for further explanations. 

quite well, although it accounts for only 2/3 of the total k. We can therefore conclude 
that the linear term (e + k)T exists in RMF whenever the isovector potential is present. 
The self-consistency between the isovector and isoscalar potentials roughly contributes to 
another 1/3. 

As shown in Fig. 3, the nuclear symmetry energy for finite nuclei calculated in RMF theory 
with PK1 is in good agreement with the experiment. The symmetry energy coefficient, as 
obtained from finite nuclei can be extrapolated to infinite nuclear matter by setting the 
surface term of Eq. (2) or Eq. (9) to zero. The theoretical asymptotic value obtained in this 
manner a' sym = A*a sym /A ps 33.6 MeV is rather close to the empirically accepted values but 
considerably smaller [by ~10%] than the so called infinite symmetric nuclear matter(NM) 
values which are equal to 37.4 MeV for NL3 [18], 36.9 MeV for TM1 [19], and 37.6 MeV for 
PK1 [20], respectively. The infinite NM values are obtained using the following formula [27]: 



sym 



2 \ m. 



rn 



*2 



(12) 



where po is t ne saturation point, kp = (37r 2 p /2) 1,/3 is the Fermi momentum and m* is the 
effective mass at saturation point. One should note that the symmetry energy coefficient a' sym 
decreases with increasing N/Z ratio in infinite NM. Moreover, there are effective Lagrangians 
having smaller values of a' sym in the infinite symmetric NM including GL-97 (32.5 MeV) [27], 
TW-99 (32.77 MeV) [28] and DD-ME1 (33.06 MeV) [29]. At present, the origin of this 
discrepancy is not clear to us. Certainly, it reflects the important role played by the nuclear 
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surface in finite nuclei, where one can notice the large contribution to the surface energy 
coming from the isovector potential. Differences in the linear term affect the size of the 
surface energy coefficient obtained in the calculations. However, the volume term of the 
total symmetry energy coefficient is not affected much even when varying the linear term 
between 1 and 4. 

In summary, the nuclear symmetry energy has been studied in RMF theory with effective 
Lagrangians NL3, TM1, and PK1. The mean level spacing e, the effective isovector potential 
strength k, and the nuclear symmetry energy coefficient a sym are calculated for the isobaric 
chains A=40, 48, 56, 88, 100, 120, 140, 160, 164, and 180 from T = to the vicinity of the 
drip line. It is shown that, except some strong variations at small values of T z , the mean 
level spacing e is stabilized at large T z and lies in the region of the empirical value after 
being re-scaled by the effective mass. These results confirm the general formulation of the 
symmetry energy obtained from the simple iso-cranking model as first proposed in Ref. [1] 
and are in agreement with Skyrme-Hartree-Fock calculations [11, 13]. 

By switching on the isovector potential due to the p meson field, the effective isovector 
potential strength k is extracted. It is surprising to find that the RMF theory, which is a 
Hartree approximation at first glance, generates a large linear term corresponding at least to 
E T — E T ~ ]-k(A)T(T+ 1 + e/k). This is in contrast to the SHF model, where the isovector 

potential has a -k(A)T(T + 1) dependence. Hence, the total nuclear symmetry energy in 
RMF follows the T(T + 1) relation quite well. 

The nuclear symmetry energy coefficient a sym as extracted from E sym = a sym T(T + 1) 
is in good agreement with the empirical data in Refs. [8, 10]. The discrepancy between 
the calculated asymptotic value of a sym and the infinite symmetric nuclear matter estimate 
requires further systematic studies. Our work also indicates that a general formulation of the 
nuclear binding energy in terms of the density functional theory in fact can yield a T(T + 1) 
dependence of the symmetry energy. The question of the physical mechanism leading to the 
restoration of the complete linear term within RMF theory is left open. 
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